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0\ . ABSTRACT 

A model for quantized gravitation based on the simplicial lattice discretization is studied in detail using 

r- 1 | 

a comprehensive finite size scaling analysis combined with renormalization group methods. The results are 
consistent with a value for the universal critical exponent for gravitation v = 1/3, and suggest a simple 

r> . 

relationship between Newton's constant, the gravitational correlation length and the observable average 

\ 

space-time curvature. Some perhaps testable phenomenological implications of these results are discussed. 
To achieve a high numerical accuracy in the evaluation of the lattice path integral a dedicated parallel 
machine was assembled. 
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1 Introduction 



One of the outstanding problems in theoretical physics is a determination of the quantum-mechanical 
properties of Einstein's rclativistic theory of Gravitation. Approaches based on linearized perturbation 
methods have had moderate success so far, as the underlying theory is known not to be perturbatively re- 
normalizable jl], ^ . Due to the complexity of even such approximate calculations a fundamental coupling of 
the theory, the bare cosmological constant term, is usually set to zero thus further restricting the potential 
physical relevance of the results. In addition gravitational fields are themselves the source for gravitation 
already at the classical level, which leads to the problem of a intrinsically non-linear theory, where pertur- 
bative results are possibly of doubtful validity for sufficiently strong effective couplings. This is especially 
true in the quantum domain, where large fluctuations in the gravitational field appear at short distances. In 
general nonpcrturbative effects can give rise to novel behavior in a quantum field theory, and in particular 
to the emergence of non-trivial fixed points of the renormalization group (a phase transition in statistical 
mechanics language). It has been realized for some time that in general the universal low and high energy 
behavior of field theories is almost completely determined by the fixed point structure of the renormalization 
group trajectories ||. 

The situation described above bears some resemblance to the theory of strong interactions, Quantum 
Chromodynamics. Non-linear effects are known here to play an important role, and end up restricting the 
validity of perturbative calculations to the high energy, short distance regime, where the effective gauge 
coupling can be considered weak due to asymptotic freedom Q . For low energy properties Wilson's discrete 
lattice formulation, combined with the renormalization group and computer simulations, has provided so far 
the only convincing evidence for quark confinement and chiral symmetry breaking, two phenomena which 
are invisible to any order in the weak coupling, perturbative expansion. 

A discrete lattice formulation can be applied to the problem of quantizing gravitation. Instead of con- 
tinuous metric fields, one deals with gravitational degrees of freedom which live only on discrete space-time 
points and interact locally with each other. In Regge's simplicial formulation of gravity ^ one approximates 
the functional integration over continuous metrics by a discretized sum over piecewise linear simplicial ge- 
ometries H 0, ||, ^1 . In such a model, the role of the continuum metric is played by the edge lengths of the 
simplices, while curvature is described by a set of deficit angles, which can be computed via known formulae 
as functions of the given edge lengths. The simplicial lattice formulation of gravity is locally gauge invariant 
and can be shown to contain perturbative gravitons in the lattice weak field expansion || , making it an 
attractive and faithful lattice regularization of the continuum theory. 
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The discretized theory is restricted to a finite set of dynemical variables, once a set of suitable boundary 
conditions are imposed such as periodic or with some assigned boundary manifold. In the end the original 
continuum theory of gravity is to be recovered as the space-time volume is made large and the fundamental 
lattice spacing of the discrete theory is sent to zero. Possibly without having to rely, at least in principle, 
on any further approximation to the original continuum theory. 

Quantum fluctuations in the underlying geometry are represented in the discrete theory by fluctuations 
in the edge lengths, which can be modeled by a well-defined, and numerically exact, stochastic process. 
In analogy with other field theory models studied by computer, calculations are usually performed in the 
Euclidean imaginary time framework, which is the only formulation amenable to a controlled numerical study, 
at least for the immediate foreseeable future. The Monte-Carlo method, based on the concept of importance 
sampling, is well suited for evaluating the discrete path integral for gravity and for computing the required 
averages and correlation functions. By a careful and systematic analysis of the lattice results, the critical 
exponents can be extracted, and the scaling properties of invariant correlation functions determined from 
first principles. 

Studies on small lattices suggest a rich scenario for the ground state of quantum gravity j?], ||, [ll| |l2| . The 
present evidence indicates that simplicial quantum gravity in four dimensions exhibits a phase transition (in 
the bare coupling G) between two phases: a strong coupling phase, in which the geometry is smooth at large 
scales and quantum fluctuations in the gravitational field eventually average out and are bounded; and a 
weak coupling phase, in which the geometry is degenerate and space-time collapses into a lower-dimensional 
manifold, bearing some physical resemblance to a branched polymer. Only the smooth, small negative 
curvature and thus anti-DeSitter-like phase appears to be physically acceptable. Phrased in different terms, 
the two phases of quantized gravity found in ]l2| , can loosely be described as having in one phase (with 
bare coupling G < G c , the rough branched polymer- like phase) 

(9n») = , (1.1) 

while in the other (with bare coupling G > G Cl the smooth phase), 

(g» v ) » C77 M „ , (1.2) 

with a vanishingly small negative average curvature in the vicinity of the critical point at G c . The existence 
of a phase transition at finite coupling G, usually associated in quantum field theory with the appearance 
of an ultraviolet fixed point of the renormalization group, implies in principle non-trivial, calculable non- 
perturbative scaling properties for correlations and effective coupling constants, and in particular in the case 
at hand for Newton's gravitational constant. Since only the smooth phase with G > G c has acceptable 
physical properties, one would conclude on the basis of fairly general renormalization group arguments that 
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at least in this lattice model the gravitational coupling can only increase with distance. Furthermore, the 
rise of the gravitational coupling in the infrared region rules out the applicability of perturbation theory to 
the low energy domain, to the same extent that such an approach is deemed to be inapplicable to study the 
low-energy properties of asymptotically free gauge theories. 

It is a remarkable property of quantum field theories that a wide variety of physical properties can 
be determined from a relatively small set of universal quantities [ fl3| . Namely the universal leading critical 
exponents, computed in the vicinity of some fixed point (or fixed line) of the renormalization group equations. 
In the lattice theory the presence of a fixed point or phase transition is often inferred from the appearance 
of non-analytic terms in invariant local averages, such as for example the average curvature 

< I 2 > < Sf x f R ( x l > ^ n(k) ~ -A n (K - k)^ , (1.3) 
< J era; y/g > fc^fcc 

where k = l/8nG. From such averages one can determine the value for is, the correlation length exponent, 

£(fc) ~ A* (fc c - k)~ v . (1.4) 

k— >k c 

An equivalent result, relating the quantum expectation value of the curvature to the physical correlation 
length £ , is 

~ £ 1/U ~ A ■ (1-5) 
Matching of dimensionalities in these equations is restored by supplying appropriate powers of the ultraviolet 
cutoff, the Planck length lp = \[G, The exponent v is known to be related to the derivative of the beta 
function for G in the vicinity of the ultraviolet fixed point, 

0'(G C ) = -l/v . (1.6) 

In addition, the correlation length £ itself determines the long-distance decay of the connected, invariant 
two-point correlations at fixed geodesic distance d. For the curvature correlation one has for distances much 
larger compared to the correlation length 

< R{x) R(y) S(\x -y\-d)> c ~ d-° e- d /Z , (1.7) 

a > ? 

while for shorter distances one expects a slower power law decay 

<^gR(x)^gR(y)6(\x-y\-d)> Cd ~^ d2{ tiM ■ ( L8 ) 

The possibility of non-trivial scaling dimensions in the theory of gravitation is not new and was pointed out 
some time ago in a series of interesting papers [|4j . Moreover it is easy to see that the scale dependence of 
the effective Newton constant is given by 

G(r) = G(0) [l + c{r/O l/v + 0{{r/d) 2/v )] , (1-9) 
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with c a calculable numerical constant. In this last expression the momentum scale £ _1 plays a role similar 
to the scaling violation parameter A-p^ of QCD. It seems natural, although paradoxical at first, to associate 
£ with some macroscopic cosmological length scale, such as the Hubble distance cHq 1 , with the lack of 
screening of gravitational interactions ultimately accounting for such an unusual interpretation Jl2| , |l5{ . Of 
course an increase of the gravitational coupling at large distances signals a likely breakdown of perturbation 
theory for computing low energy properties of gravity. 

It should be clear, even from this brief discussion, that the critical exponents by themselves already 
provide a significant amount of useful information about the continuum theory. In reality, the complexity of 
the lattice interactions and the practical need to sample many statistically independent field configurations 
contributing to the path integral, which is necessary for correctly incorporating into the model the effects 
of quantum-mechanical fluctuations, leads to the requirement of powerful computational resources. The 
results presented in this paper were obtained using a dedicated custom-built 20-GFlop 64-processor parallel 
computer, described in detail in p6| . 

Finally one should mention that recently there has been a significant resurgence of interest in the classical 
applications of the Regge formulation to gravity. A description of the methods as applied to several aspects 
of the initial value problem in General Relativity can be found in the recent references in |rj| . For a related 
approach to lattice gravity based on dynamical triangulations see also . 

A brief outline of the paper is as follows. Section 2. contains a discussion of general and finite size scaling 
and related issues as they apply to the lattice theory of gravity. Section 3. touches on the issue of the 
unboundedness of the Euclidean gravitational action. Section 4. defines local curvature averages and their 
fluctuations, while Section 5. introduces a set of exact sum rules for averages which follow from the scaling 
properties of the partition function. Section 6. defines a set of invariant correlations and discusses how they 
relate to the local fluctuations defined previously. Section 7. includes a general discussion of the expected 
properties of the theory in the presence of an ultraviolet fixed point, including expectations based on the 
analytical 2 + e expansion. In Section 8. the numerical results are presented. Section 9. contains a discussion 
of the possible future physical relevance of the results, while Section 10. contains the conclusions. 



2 Finite Size Scaling 



One of the most important quantities used in establishing the continuum limit of a lattice field theory are 
the critical exponents. Reliable estimates for the exponents in a lattice field theory require a comprehensive 
finite-size analysis, a procedure by which accurate values for the critical exponents are obtained by taking 
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into account the linear size dependence of the result computed in a finite volume V. One starts from the 
general Euclidean Action (or statistical mechanics Hamiltonian) 

H = 9* Oi (2.10) 

i 

with <?j the coupling associated with the operator Oi. In the gravitational case the couplings would correspond 
to the bare cosmological constant, the Newtonian gravitational constant and the higher derivative coupling. 
Close to a renormalization group fixed point denoted by {g*} one chooses the Oi's to be eigenvectors of the 
linearized renormalization group transformation, such that 

9i-9t "> Wifr-gt) , (2.11) 

where b is the scale factor of the transformation. In the simplest statistical mechanics systems, such as a 
fcrromagnet in the absence of an external magnetic field, one has O ~ H as the only relevant operator (in 
the sense that y > 0), and g ~ t = T — T c . As will be discussed below, in the gravitational case the role of 
T is played by the bare gravitational coupling G. Additional operators appearing in the action are classified 
as marginal (y = 0) or irrelevant. The relevance of the energy operator reflects the fact that close to the 
critical point t is the only parameter that needs to be tuned to achieve criticality, synonymous with long 
range correlations. Universality of critical behavior then accounts for the fact that many diverse physical 
systems exhibit the same scaling behavior in the vicinity of the critical point, as a consequence of a divergent 
correlation length |H| . 

In practice the renormalization group approach is brought in via a slightly different route, involving a 
change in the overall linear size of the system. The usual starting point for the derivation of the scaling 
properties of the theory is the Renormalization Group (RG) behavior of the free energy F = — log Z/V 

F(t,{ Uj }) = F reg (t,{u 3 }) + b- d F sing (b y %{b^ Uj }) , (2.12) 

where F S i ng is the singular, non-analytic part of the free energy, and F reg is the regular part, b is the block 
size in the RG transformation, while y t and y 3 \j > 2) are the relevant eigenvalues of the RG transformation 
(for more details see the review |fj^]). One denotes here by y t > the relevant eigenvalue, while the remaining 
eigenvalues yj < are associated with either marginal or irrelevant operators. Usually y t _1 is called v, while 
the next subleading exponent y<i is denoted —lj. 

The correlation length £ determines the asymptotic decay of correlations, in the sense that one expects 
for example for the two-point function at large distances 

< 0{x)0(y) > - e -\*-v\/t . (2.13) 

\x-y\ > £ 
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The scaling equation for the correlation length itself 

C(t) = (2.14) 

implies for & = t~ l ' Vt that £ ~ t~ v with a correlation length exponent 

v = l/yt ■ (2.15) 

Derivatives of the free energy F with respect to t then determine, after setting the scale factor b — i _1 / yt , 
the scaling properties of physical observables, including corrections to scaling (2t|. Thus for example, the 
second derivative of the free energy with respect to t yields the specific heat exponent a = 2 — d/yt = 2 — dp, 

^F(t,{u 3 }) ~ f -0-*O . (2.16) 

In the gravitational case one identifies the scaling field i with fc c — k, where k = l/16nG involves the 
bare Newton's constant. The appearance of singularities in physical averages, obtained from appropriate 
derivatives of F, is rooted in the fact that close to the critical point at t = the correlation length diverges. 

The above results can be extended to the case of a finite lattice of volume V and linear dimension 
L = V x l d . The volume-dependent free energy is then written as 

Ffafa^L- 1 ) = F reg {t,{u 3 }) + b- d F sing (bv%{bW Uj },b/L) . (2.17) 

For b = L (a lattice consisting of only one point) one obtains the Finite Size Scaling (FSS) form of the free 
energy (for a detailed presentation of this procedure see ; see also |^2|, for a field-theoretic justification) . 
After taking derivatives with respect to the fields t and {%}, the FSS scaling form for physical observables 
follows. For a quantity O diverging like t~ xo in the infinite volume limit one has 



0(L,t) = L 



_ TX /V 



(2.1S) 



with fo a smooth scaling function, and £(oo,i) the infinite volume correlation length. For sufficiently large 
volumes the correction to scaling term involving lo can be neglected, but in general one needs to be aware 
of their presence if either the volumes are not large enough or if the corrections are large due to a large 
amplitude or small exponent. Some properties of the scaling function fo{y) can be deduced on general 
grounds: it is expected to show a peak if the finite volume value for O is peaked, it is analytic at x = since 
no singularity can develop in a finite volume, and fo(y) ~ V ~ x ° for large y for a quantity O which diverges 
as t~ x ° in the infinite volume limit. 

The last expression is useful when the infinite-volume correlation length is known. But since close to the 
critical point £ ~ t~" ' , one can deduce the equivalent scaling from 

0(L,t) - L*°' v \fo(Lf) + 0(L- u j\ , (2.19) 
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which relies on a knowledge of i, and thus of the critical point, instead. For a state of the art application of 
the above methods to the 3-d Ising model see 

The previous discussion applies to continuous, second order phase transitions. First order phase transi- 
tions are driven by instabilities, and are in general not governed by any renormalization group fixed point. 
The underlying reason is that the correlation length does not diverge at the transition point, and thus the 
system never becomes scale invariant. Exponents for continuous, second order phase transition in general 
obey the rigorous bound 

y t < d or v > 1/d . (2.20) 

First order phase transition in renormalization group theory, on the other hand, can be associated with the 
somewhat pathological case v = 1 /d, for which the first derivative of the free energy develops a step-function 
singularity. In a renormalization group framework the corresponding pseudo-critical point is denoted as a 
discontinuity fixed point |25| . 

In the simplest case, a first order transition develops as the system tunnels between two neighboring 
minima of the free energy. In the metastable branch the free energy acquires a complex part with an 
essential singularity in the coupling located at the first order transition point [^6[ ^7|. As a consequence, such 
a singularity is not generally visible from the stable branch, in the sense that a power series expansion in the 
temperature is unaffected by such a singularity. Indeed in the infinite volume limit the singularity associated 
with a first order transition at T c becomes infinitely sharp, like a S— or ^-function type singularity. The 
singularity in the free energy at the endpoint of the metastable branch (at say T*) then cannot be explored 
directly, it has to reached by an analytic continuation from the stable side of the free energy branch. 



3 Unboundedness of the Euclidean Theory 



Perturbation theory on a lattice and in the continuum suggests the presence of an instability in the Eu- 
clidean formulation for sufficiently smooth manifold. It is also known that the above instability is associated 
with the appearance of a wrong sign for the conformal mode. On the lattice the instability seems to persist 
close to the critical point (l2), which suggests that the continuum limit has to be reached by some sort 
of analytic continuation from the stable phase towards the critical point, naturally defined as the point in 
coupling constant space where the correlation length diverges. 

In the weak-field expansion psj the Einstein-Hilbert action contains both spin two (graviton) and spin 
zero (conformal mode) contributions. In the continuum one can by a judicious choice of invariant correlation 
functions isolate physical properties of the graviton from the conformal mode. A similar result holds on 
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the lattice, as can be seen by expanding the Regge action about a regular lattice and using the fact that 
the lattice and continuum actions are equivalent for sufficiently smooth manifolds ^9). In general, after 
expanding the metric around flat space (which requires A = 0), 



V16ttG V , (3.21) 
one can cast the lowest order quadratic contribution to the action in the form 

/e[/i M "] = \ \ <^ x h^V^xuhxa , (3.22) 



where V is a matrix which can be expressed in terms of spin projection operators. In momentum space it 
can be written as 

V = [P {2) - 2P(°)] p 2 , (3.23) 

where P^ and p(°) are spin two and spin zero projection operators introduced in poj . Physically, the two 
terms correspond to the propagation of the graviton and of the conformal mode, respectively, with the latter 
one appearing with the 'wrong' sign. In the 'Landau' gauge, with a gauge fixing term a~ l (d^h^ u ) 2 and 
a = 0, one obtains for the graviton propagator in momentum space 

p(2) lp(0) 

<W(P) = -Jf^l- . (3.24) 

The unboundedness of the Euclidean gravitational action shows up clearly in the weak field expansion, with 
the spin zero mode acquiring a propagator term with the wrong sign, i It has been argued that in weak field 
perturbation theory and in order to avoid the unboundedness problem one should perform the functional 
integral over metrics by distorting the integration contour so as to include complex conformal factors 
One drawback of this prescription is that it only appears applicable within the framework of perturbation 
theory. For a recent review of the Euclidean instability problem see |32| . 

In the presence of a cosmological constant, things are further complicated by the fact that since flat space 
is no longer a solution of the classical equations of motion, and the above expansion for the metric looses part 
of its meaning due to the presence of the tadpole term. But after shifting to the correct 0-th order solution, 
a similar result is obtained. One can further modify the action to include additional invariant terms, but 
things do not get any better. In the presence of higher derivative terms in the gravitational action, the above 
result is modified by terms 0(p 4 ), and becomes |33| 

p(2) lp(0) 

•Wri-^fS^ + ^r^ ■ (3 - 25) 



2 It should be noted that such an instability is not peculiar to gravitation. Indeed the Euclidean path integral for the one- 
dimensional Coulomb potential, an otherwise completely well behaved quantum mechanical system, already exhibits such an 
instability. It would be premature to conclude from such a result that the problem is physically ill-posed. 
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The p 4 terms improve the ultraviolet behavior of the theory, but do not remove the unboundedness problem, 
which re-appears for sufficiently small p 2 , in the low momentum or long-distance limit. Moreover, the 
resulting theory is most likely not unitary unless the coupling a is vanishingly small. The lack of positivity 
of physical correlations for a > can be seen explicitly even in a non-perturbative treatment pi| , and makes 
such a modified theory of gravitation in the end somewhat unattractive. 



4 Local Averages and Fluctuations 



In the following the relevant definitions for gravitational averages and correlations on the lattice will be 
briefly recalled, in a form which will be used in later sections. The starting point for a non-perturbative 
study of quantum gravity is a suitable definition of the discrete Feynman path integral. In the simplicial 
lattice approach one starts from the discretized Euclidean path integral for pure gravity, with the squared 
edge lengths taken as fundamental variables, 

Zl = f]] (W)* II dl l ex P | - E ( A Vh - k ShAh + a % A h/Vh +■■■)) • (4-26) 

"' s ij { h J 

The above expression represents a lattice discretization of the continuum Euclidean path integral for pure 
quantum gravity 
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j n (V^) CT n d ^ x ) exp {-/ d " x ^ i x - ^R+i R ^p* Rfivpc +■••)} > ( 4 - 27 ) 

with k^ 1 — 8ttG, and G Newton's constant, and reduces to it for smooth enough field configurations. In 
the discrete case the integration over metrics is replaced by integrals over the elementary lattice degrees 
of freedom, the squared edge lengths. The discrete gravitational measure in Zl can be considered as the 
lattice analog of the DeWitt |55| continuum functional measure Jl5|] . The 8 A term in the lattice action is 
the well-known Regge term and reduces to the Einstein-Hilbert action in the lattice continuum limit 
H p9[ . A cosmological constant term is needed for convergence of the path integral, while the curvature 
squared term allows one to control the fluctuations in the curvature @, ^, |ll, 12 . In practice, and for obvious 
phcnomenological reasons, one is only interested in the limit when the higher derivative contributions are 
small compared to the rest, a — > 0. In this limit the theory depends, in the absence of matter and after a 
suitable rescaling of the metric, only on one bare parameter, the dimcnsionlcss coupling k 2 /X. Without loss 
of generality, one can therefore set the bare cosmological constant A = 1 . 

Some partial information about the behavior of physical correlations can be obtained indirectly from 
local invariant averages. In (?], [l2| gravitational observables such as the average curvature and its fluctuation 
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were introduced. The appropriate lattice analogs of these quantities are readily written down by making use 
of the usual correspondences J d A x ^fg — > Xhingcsh Vh etc.. On the lattice the natural choices for invariant 
operators are 

hinges hZ)x 

^gR(x) -> 2 ]T 5 h A h 

hinges h~2>x 

y/gR^XaBT^ix) -> 4 ^ ((SfcAfc) 2 /^ (4.28) 

hinges ^Dx 

(we have omitted here on the r.h.s. an overall numeric coefficient, which will depend on how many hinges 
are actually included in the summation; if the sum extends over all hinges within a single hypercube, then 
there will be a total of 50 hinge contributions). In this paper no higher derivative terms will be considered, 
and thus only the first and second operators will be used in the following discussion. 

On the lattice one prefers to define quantities in such a way that variations in the average lattice spacing 



V < I 2 > are compensated by the appropriate factor as determined from dimensional considerations. In the 
case of the average curvature one defines the lattice quantity 1Z as 



K(k) ee < I* > <2E.^> j (4 29) 



which in the continuum corresponds to 

> • /.• i 

< / d 4 x ^fg > 



_ < Sf '^m > . (430) 



and similarly for the curvature fluctuation, 



XR(k) _ <E>W)- . (431) 

<2^h V h> 



which in the continuum corresponds to 



„<», „ <U*s>>>-<^*>' . (432) 



The latter is related to the connected curvature correlation at zero momentum 



fd*xfd*y < ^gjx)R(x) y/^y)R(y) > c 

Xn f— — r= • {l.SS) 

< J d 4 xy/g{x) > 



Both 1Z and \n are related to derivatives of Zl with respect to k, 



and 



n(k) ~ V§k lnZL ' (434) 



Xn(k) ~ (4.35) 
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One can contrast the behavior of the preceding quantities, associated strictly with the curvature, with the 
analogous quantities involving the local volumes (and which correspond to the square root of the determinant 
of the metric in the continuum). Consider the average volume per site 



and its fluctuation defined as 



No 

h 



xm . <(S>^-<E^> 2 . (437) 

where one denotes by Vh the volume associated with the hinge h. In the continuum it corresponds to the 
expression 

%m „ <'/^>J_< > /v^> 2 . (438) 

The latter is related to the connected volume correlator at zero momentum 

^ JjMjWlyfe (4 39) 
< / d i x v r ^(x) > 

The average volume per site < V > and its fluctuation xv are simply related to derivatives of Zl with respect 
to the bare cosmological constant A, 

<V>~ j^lnZ L , (4.40) 

and 

Xv(k) ~ ^\nZ L . (4.41) 

One would expect the fluctuations in the curvature to be sensitive to the presence of a spin two massless 
particle, while fluctuations in the volume would only probe the correlations in the conformal mode channel. 



5 Sum Rules 



In this section some useful sum rules will be derived, which follow from simple scaling properties of the 
discrete functional integral. These will be later used in the discussion of the numerical results. A simple 
scaling argument, based on neglecting the effects of curvature terms entirely (which vanish in the vicinity of 
the critical point), gives first of all an estimate of the average volume per edge 

<^>~ W+f* ~ * (5.42) 

\d d=4,a=0 2A V ^ 

In four dimensions the numerical simulations with a = agree quite well with the above formula. 
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Additional exact lattice identities can be obtained by examining the scaling properties of the action and 
measure. The bare couplings k and A in the gravitational action are dimcnsionful in four dimensions, but 
one can define the dimensionless ratio k 2 /A, and rescale the edge lengths so as to eliminate the overall length 
scale yJk/X. As a consequence the path integral for pure gravity, 



Z L {X,k,a) = / dfi[l 2 } e- I[l2] 



obeys the simple scaling property 



' k \ Nl /k 2 k 2 

Z L (\.k.,n= [j) ^ ppa] = (A, 



-Ni/2 



z, [ i,^=,a) = (xy Nl z L (A l,a 



(5.43) 



(5.44) 



where Ni represents the number of edges in the lattice, and the dl 2 measure (a = 0) has been selected fl5| , 
which is the lattice analog of the continuum DeWitt functional measure. This equation implies in turn a 
sum rule for local averages, which (again for the specific case of the dl 2 measure) reads 



2A < J2 V h > - k <J2 ShAh > - Ni = , 



(5.45) 



and is easily derived from Eq. ( 5.44 ) and the definitions in Eqs. ( 4.34 ) and ( 4.40| ). No represents the number 
of sites in the lattice, and the averages are defined per site (for the hypercubic lattice used in this paper, 
Ni = 15A , N 2 = 50Ao, N 3 = 36N and N A = 24iV ). The coefficients on the l.h.s. of the equation reflect 
the scaling dimensions of the various terms, with the last term on the l.h.s. term arising from the scaling 
property of the functional measure. This last formula is very useful in checking the accuracy of numerical 
calculations and the convergence properties of the Monte Carlo sampling, and is usually satisfied to high 
accuracy O(10 -4 ). It is easy to see that a similar sum rule holds for the fluctuations, 



AX 2 



< 



2 > - < J2 Vh > 2 



< 



2 >-<Y. ShAh >2 



2Ni = 



(5.46) 



h h 

Further sum rules can be derived by considering even higher derivatives of lnZj, with respect of A and k. 
The last equation relates the fluctuation in the curvature to fluctuations in the volumes, and thus implies 
a relationship between their singular parts as well. In particular, a divergence in the curvature fluctuation 
implies a divergence of the same nature in the volume fluctuation. In light of the previous discussion, from 
now on we shall consider without loss of generality only the case of bare coupling A = 1. As a consequence, 
all lengths will be tacitly expressed in units of the fundamental microscopic length scale A -1 / 4 . 



6 Invariant Correlations 
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In quantized gravity complications arise due to the fact that the physical distance between any two points 
x and y in a fixed background geometry, 

r(y) 



d{x,y\g) = min / drJg^^tfr > ( 6 - 47 ) 

« Jt(x) 

is a fluctuating quantity dependent on the choice of background metric. In addition, the Lorentz group used 
to classify spin states is meaningful only as a local concept. Since the simplicial formulation is completely 
coordinate independent, the introduction of the local Lorentz group requires the definition of a tetrad within 
each simplex, and the notion of a spin connection to describe the parallel transport of tensors between flat 
simplices. Some of these aspects have recently been discussed from a continuum point of view in J|7], [38| |39fl . 

If the deficit angles are averaged over a number of contiguous hinges which share a common vertex, one 
is naturally lead to the connected correlator 

G R (d) =<J2 § h A h Sh ' Ah ' S ^ x ~y\~ d )>c (6.48) 

h~Dx h'Dy 

which probes correlations in the scalar curvatures 

G R (d) ~ < V5 R(x) R(y) S(\x ~y\~d)> c . (6.49) 
Similarly one can construct the connected correlator 

G v (d) = < V h J2 Vh ' S (\ x -v\-d)>c, (6.50) 

hDx h'Dy 

which probes correlations in the volume elements 

G v - < y/g{x) V?(y) «5(|a: - y\ - d) > c . (6.51) 

The correlation length £ is defined through the long-distance decay of the connected, invariant correlations 
at fixed geodesic distance d. For the curvature correlation one has, at large distances, 

< y/gR(x) V9R(y)8(\x-y\-d)> c ~ . (6.52) 

d > 4 



i \ 2n 



At shorter distances one expects a slower, power law decay 

<^R(x)^R(y)S(\x-y\-d)> c I ~ j (6.53) 

with a power characterized by the exponent n. In both cases, the distances considered are much larger than 
the lattice spacing, d, £ ^> Iq. From scaling considerations one can show (see below) n = 4— 1/v. 

Simple scaling arguments allow one to determine the scaling behavior of correlation functions from the 
critical exponents which characterize the singular behavior of local averages in the vicinity of the critical 
point. A divergence of the correlation length £ 

f(Jfe) = m(k)- 1 ~ A^{k c -k)- v (6.54) 
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signals the presence of a phase transition, and leads to the appearance of a singularity in the free energy 
F(k). The presence of a phase transition usually inferred from non-analytic terms in invariant averages, such 
as the average curvature. The curvature critical exponent S is introduced via 

TZ(k) ~ -A n (k c -k) s . (6.55) 

k — >k c 

An additive constant could be added, but the evidence up to now points to this constant being zero. Similarly 
one sets for the curvature fluctuation 

Xn{k) ~ ~A n (k c ~k)- (1 -V . (6.56) 

Scaling (Eqs. 2.15| ) relates the exponent 6 to v, 

v = i±* . (6.57) 
a 

From such averages one can determine the value for v, the correlation length exponent, An equivalent result, 
relating the quantum expectation value of the curvature to the physical correlation length £ , is 

K(Z) t ~ e /u -* > (6-58) 

k — >k c 



which is obtained from Eqs. ( 6.54 ) and fl6.55 ) using 6.57 ). Matching of dimensionalities in these equations 



is restored by supplying appropriate powers of the Planck length lp — VG. 

It is then easy to relate the critical exponent v to the scaling behavior of correlations at large distances. 
The curvature fluctuation is related to the connected scalar curvature correlator at zero momentum 

. MjilMM^ (h _ „ r * . (6 . 59) 

< J d^x^fg > k-*k c 
A divergence in the fluctuation is then indicative of long range correlations, corresponding to the presence 
of a massless particle. Very close to the critical point one would expect for large separations a power law 
decay in the geodesic distance, 



1 

\x-y\->oo [x-y^ 



< ^R(x)^R(y) > t _ _~ u , (6.60) 



with the power n related to the exponent 5 via n = Sd/(1 + 5) = d — \ jv. A priori one cannot exclude to 
possibility that some states acquire a mass away from the critical point, in which case one would expect the 
following behavior for the correlation functions, 

<yfgR{x)y/gR( y )> ~ exp(-|x-y|/0 , (6.61) 

where £ is the fundamental correlation length, and m — l/£ the associated mass. The above equation can 
in fact be considered as a definition for what is meant by the correlation length £. 
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7 Beta function and Continuum Limit 



The long distance behavior of quantum field theories is determined by scaling behavior of the coupling 
constant under a change in the momentum scale. Asymptotically free theories such as QCD lead to vanishing 
gauge couplings at short distances, while the opposite is true for QED. In general the fixed point(s) of the 
renormalization group need not be at zero coupling, but can be located at some finite G c , leading to a 
non-trivial fixed point or limit cycle (|[ [IJ j40| . 

In the 2 + e perturbative expansion for gravity fi"l|| one analytically continues in the spacetime dimension 
by using dimensional regularization, and applies perturbation theory about d = 2, where Newton's constant 
is dimensionless. A similar method is quite successful in determining the critical properties of the 0(n)- 
symmetric non- linear sigma model above two dimensions p2| . In this expansion the dimensionful bare 
coupling is written as Go = A 2 ~ d G, where A is an ultraviolet cutoff (corresponding on the lattice to a 
momentum cutoff of the order of the inverse average lattice spacing, A ~ l//o)' There seem to be some 
technical difficulties with this expansion due to the presence of kinematic singularities for the graviton 
propagator in two dimension (the Einstein action is a topological invariant in d = 2), but which seem to 
have been overcome recently. A double expansion in G and e = d— 2 then leads in lowest order to a nontrivial 
fixed point in G above two dimensions 

P(G) = = (d-2)G- /3 G 2 +••• , (7.62) 

with Po > for pure gravity. To lowest order the ultraviolet fixed point is then at G c — l/(3a(d — 2). 



Integrating Eq. ( 7.62 ) close to the non-trivial fixed point one obtains for G > G c 



m = K ™[-fwk)°^ AlG - G < rmo ° } ■ (7 - 63) 

where to is an arbitrary integration constant, with the dimensions of a mass, and which should be associated 
with some physical scale. It would appear natural here to identify it with the inverse of the gravitational 
correlation length (£ = to -1 ), or some scale associated with the average curvature. The derivative of the 
beta function at the fixed point defines the critical exponent v, which to this order is independent of /3q, 
(3'(G C ) = -(d-2) = -1/v. 

The previous results illustrate how the lattice continuum limit should be taken. It corresponds to A — > oo, 
G — ► G c with to held constant; for fixed lattice cutoff the continuum limit is approached by tuning G to G c . 
In four dimensions the exponent v is defined by 

m ~ Ck\G-G c \ v , (7.64) 



16 



where m is proportional to the graviton mass, and C is a calculable numerical coefficient. The value of 
v determines the running of the effective coupling G(/i), where is an arbitrary momentum scale. The 
renormalization group tells us that in general the effective coupling will grow or decrease with length scale 
?' = l//i, depending on whether G > G c or G < G Cl respectively. The physical mass parameter m is itself 
scale independent, and obeys the Callan-Symanzik renormalization group equation 

»Q^rn = n— {Cii\G{ii)-G c \ v } = . (7.65) 

As a consequence, for G > G c , corresponding to the smooth phase, one expects for the running, effective 
gravitational coupling Jl^, |1| 

G(r) = G(0) [ 1 + c(r/0 1/u + O((r/0 2/u ) ] , (7.66) 

with c a calculable numerical constant. The physical mass m = £ determines the magnitude of scaling 
corrections, and plays a role similar to A-p^ in QCD. ft cannot be determined perturbatively as it appears 
as an integration constant. Physically it separates the short distance, ultraviolet regime with characteristic 
momentum scale fi, 

l^ 1 > jti > m , (7.67) 

from the large distance, infrared region 

m>/j>i" 1 , (7.68) 

where L =<V> 1 ^ 4 is the linear size of the system. 

The exponent v is simply related to the derivative of the beta function for G in the vicinity of the 
ultraviolet fixed point, 

/3'(G C ) = -1/p . (7.69) 

Thus computing v is equivalent to computing the derivative of the beta function in the vicinity of the 
ultraviolet fixed point. There are indications from the lattice theory that only the smooth phase with 
G > G c exists (in the sense that spacetime collapses onto itself for G < G c ), which would suggest that the 
gravitational coupling can only increase with distance. 

One should also perhaps recall here the fact that a bare cosmological constant A, which could appear in 
the original action (as indicated in Eq. ( 4.26| )) has been scaled out, when it was set equal to one by rescaling 



all the edge lengths. If one puts it back in, then the effective Newton's constant would have to be multiplied 
by that bare scale. As a result one obtains for the running of Newton's constant, valid for "short" distances 
(i > m, 

GW r<« P ° X ~ 1/2 Gc+ lef ' _ ' ' (7 - 7()l 
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where G c is a pure number of order one, and below it will be argued that 1/v = 3. The quantity Iq is the 
average lattice spacing, and the correct dimensions for G(/i) (length squared) have been restored. In addition 



a bare cosmological constant A was re-introduced, which was previously set equal to one in Eq. (4.26) (it 
fixes the overall length scale in the functional integral over edge lengths). 



8 Numerical Results 



Next we come to a discussion of the numerical methods employed in this work and the analysis of the 
results. As in previous work, the edge lengths are updated by a straightforward Monte Carlo algorithm, 
generating eventually an ensemble of configurations distributed according to the action and measure of 
Eq. (4.26). Further details of the method as applied to pure gravity are discussed in |?], [l^], and will not be 



repeated here. 

In this work lattices of size 4x4x4x4 (with 256 sites, 3840 edges, 6144 simplices) 8x8x8x8 (with 
4096 sites, 6144 edges, 98304 simplices) 16 x 16 x 16 x 16 (with 65536 sites, 983040 edges, 1572864 simplices) 
were considered. Even though these lattices are not very large, one should keep in mind that due to the 
simplicial nature of the lattice there are many edges per hypercube with many interaction terms, and as 
a consequence the statistical fluctuations can be comparatively small, unless measurements are taken very 
close to a critical point, and at rather large separation in the case of the potential. The results presented 
here are still preliminary, and in the future it should be possible to repeat such calculations with improved 
accuracy on much larger lattices. 

The topology is restricted to a four-torus (periodic boundary conditions). We have argued before that 
one could perform similar calculations with lattices employing different boundary conditions or topology, 
but the universal infrared scaling properties of the theory should be determined only by short-distance 
renormalization effects. 

It seems reasonable that based on physical considerations one needs to impose the constraint that the 
scale of the curvature be much smaller than the average lattice spacing, but still much larger than the overall 
size of the system. In other words 

<l 2 > < <l 2 >\n\~ 1 < <V> 1/2 . (8.71) 

Or, that in momentum space the physical scales should be much smaller that the ultraviolet cutoff, but much 
larger than the infrared cutoff. An equivalent requirement is then 

L^ 1 < m < ^o 1 > ( 8 - 72 ) 
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where L is the linear size of the system, m = l/£, and Iq the lattice spacing. It should be kept in mind 
that in this model, and contrary to ordinary gauge theories on a lattice, the lattice spacing is a dynamical 
quantity. Even close to the critical point where the curvature vanishes the lattice is by no means regular, 



and the quantity Iq = v<P> only represents an "average" cutoff parameter. 



The bare cosmological constant A appearing in the gravitational action of Eq. (4.26) was fixed at 1 (since 
this coupling sets the overall length scale in the problem), and the higher derivative coupling a was set to 
(pure Regge-Einstein action). For the measure in Eq. ( |4.2g ) this choice of parameters leads to a well behaved 
ground state for k < k c sa 0.053 for a = [Jl2|, [ll). The system then resides in the 'smooth' phase, with 
a fractal dimension close to four; on the other hand for k > k c the curvature becomes very large ('rough' 
phase), and the lattice tends to collapse into degenerate configurations with very long, elongated simplices 
f| [ll], For a = we investigated 22 values of k. 

On the 16 4 lattice 36,000 consecutive configurations were generated for each value of k, and 22 different 
values for k were chosen. The results for different values of k can be considered as completely statistically 
uncorrelated, since they originated from unrelated configurations. On the smaller 8 4 lattice 100,000 consec- 
utive configurations were generated for each value of k. On the 4 4 lattice 500,000 consecutive configurations 
were generated for each value of k. To accumulate the results, the machine ran continuously for about 14 
months. 



The results obtained for the average curvature 1Z (defined in Eq. 4.29) as a function of the bare coupling 
k are shown in Fig. 1., on lattices of increasing size with 4 4 , 8 4 and 16 4 sites. Fig. 2. shows the 16 4 data by 
itself. The errors in are quite small, of the order of a tenth of a percent or less, and are therefore not visible 
in the graph. 

In fll2| it was found that as k is varied, the curvature is negative for sufficiently small k ('smooth' phase), 
and appears to go to zero continuously at some finite value k c . For k > k c the curvature becomes very 
large, and the simplices tend to collapse into degenerate configurations with very small volumes (< V> / < 
I 2 > 2 ^ 0). This 'rough' or 'collapsed' phase is the region of the usual weak field expansion (G — > 0); in the 
continuum it is characterized by the unbounded fluctuations in the conformal mode. But there appears to 
be more structure to the data. 

Accurate and reproducible curvature data can only be obtained for k below the instability point k u since, 
as already pointed out in |i~2| , for k > k u ss 0.053 an instability develops, presumably associated with the 
unbounded conformal mode. Its signature is typical of a sharp first order transition, beyond which the 
system tunnels into the rough, elongated phase which is two-dimensional in nature and has no physically 
acceptable continuum limit. The instability is caused by the appearance of one or more localized singular 
configuration, with a spike-like curvature singularity. It is not associated with any sort of coherent effect or 
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the appearance of long-range order, and remains localized around a few lattice points. In other words, the 
correlation length £ remains finite at k u . At strong coupling such singular configurations are suppressed by a 
lack of phase space due to the functional measure, which imposes non-trivial constraints due to the triangle 
inequalities and their higher dimensional analogs. In other language, the measure regulates the conformal 
instability at sufficiently strong coupling. 

It is characteristic of first order transitions that the free energy develops only a delta-function singularity 
at k u , with the metastable branch developing no non-analytic contribution at k u . Indeed it is well known 
from the theory of first order transitions that tunneling effects will lead to a purely imaginary contribution 
to the free energy, with an essential singularity for k > k u p6| . In the following we shall clearly distinguish 
the instability point k u from the true critical point k c . 

As a consequence, the non-analytic behavior of the free energy (and its derivatives which include for 
example the average curvature) has to be obtained by analytic continuation of the Euclidean theory into 
the metastable branch. This procedure, while unusual, is formally equivalent to the construction of the 
continuum theory exclusively from its strong coupling (small k, large G) expansion 

oo 

Z L {k) = ^a„fc" , (8.73) 

oo 

Jl(k) = ^6„fc" . (8.74) 

n=Q 

Given a large enough number of terms in this expansion, the nonanalytic behavior in the vicinity of the 
true critical point at k c can then be determined using differential or Pade approximants p3| , for appropriate 
combinations of thermodynamic functions which are expected to be meromorphic in the vicinity of the true 
critical point p4j . In the present case, instead of the analytic strong coupling expansion, one has at one's 
disposal a set of (in principle, arbitrarily) accurate data points to which the expected functional form can 
equally be fitted. And what is assumed is the kind of regularity which is always assumed in extrapolating 
finite series (whether convergent or asymptotic as in the case of QED or A</> 4 in d < 4 f45{| ) to the boundary 
of their radius of convergence. 

Ultimately it should be kept in mind that one is really only interested in the pseudo-Riemannian case, 
and not the Euclidean one for which an instability due to the conformal mode is to be expected. Indeed 
had such an instability not occurred one might wonder if the resulting theory still had any relationship to 
the original continuum theory. Arguments based on effective actions suggest that if the Euclidean (or more 
appropriately, Riemannian) lattice theory eventually approaches the classical continuum theory at large 
distances and in the vicinity of the critical point, then an instability in the quantum lattice theory must 
develop, since the continuum classical theory is known to be unstable. 
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In the following only data for k < k u will be considered; in fact to add a margin of safety only k < 0.051 
will be considered throughout the rest of the paper. This choice will avoid the inclusion in the fits of any 
data affected by the sharp turnover which appears, for large lattices, at k — k u ss 0.053. 

To extract the critical exponent 5, one fits the computed values for the average curvature to the form 
(see Eq. 6.55| ) 



K{k) -A n {k c -k) s . (8.75) 

It would seem unreasonable to expect that the computed values for 1Z are accurately described by this 
function even for small k. Instead the data is fitted to the above functional form for either k > 0.02 or 
k > 0.03 and the difference in the fit parameters can be used as one more measure for the error. Additionally, 
one can include a subleading correction 

TZ(k) -A n [k c ~k + B {k c -kf } S , (8.76) 

and use the results to further constraint the errors on A-r, k c and 5 = Av — 1. 
Using this set of procedures one obtains on the lattice with 4 4 sites 

k c = 0.0676(20) v = 0.343(8) (8.77) 

and on the lattice with 8 4 sites one finds 

k c = 0.0614(27) v = 0.322(16) (8.78) 

while on the lattice with 16 4 sites one finds 

k c = 0.0630(11) v = 0.330(6) . (8.79) 

These results suggest that v is very close to 1/3, and can be compared to the older low-accuracy estimate 
on an 8 4 lattice obtained in |l^] for a = 0, v — 0.33(3). 

Fig. 3. shows a graph of the average curvature TZ(k) raised to the third power. One would expect to 
get a straight line close to the critical point if the exponent for TZ(k) is exactly 1/3. The numerical data 
indeed supports this assumption, and in fact the linearity of the results close to k c is quite striking. The 
computed data is quite close to a straight line over a wide range of k values, providing further support for 
the assumption of an algebraic singularity for TZ(k) itself, with exponent close to 1/3. Using this procedure 
one finds on the 16 4 -site lattice 

k c = 0.0639(10) . (8.80) 

Since the critical exponents play such a central role in determining the existence and nature of the 
continuum limit, it appears desirable to have an independent way of estimating them, which either does not 
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k = 1/8-kG 



Fig. 1 . Average curvature 1Z as a function of k, on lattices with 4 4 (□), 8 4 (A) and 16 4 (o) sites. Statistical errors 
(~ 0(1CF 3 )) are much smaller than the size of the symbols. The thin-dotted, dotted and continuous lines represent 
best fits of the form TZ(k) = A (k c - k) s . 
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-0.01 0.01 0.02 0.03 0.04 0.05 0.06 0.07 

k = 1/8-kG 



Fig. 2 . Average curvature 1Z as a function of k, on the 16 4 (o) lattice only. Statistical errors (~ 0(1O~ 3 )) are much 
smaller than the size of the symbols. The continuous line represents a best fit of the form A (k c — k) s for k > 0.02, 
with 8 — 4v — 1 . 
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-0.01 0.01 0.02 0.03 0.04 0.05 0.06 0.07 

k = 1/8ttG 



Fig. 3 . Average curvature on the 16 4 lattice, raised to the third power. If 8 = v = 1/3, the data should fall on a 
straight line. The continuous line represents a linear fit of the form A (k c — k). The small deviation from linearity of 
the transformed data is quite striking. 
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k c = 0.0637 
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Fig. 4 . Finite size scaling behavior of the scaled curvature versus the scaled coupling. Here L = 4 for the lattice 
with 4 4 sites (□), L = 8 for a lattice with 8 4 sites (A), and L = 16 for the lattice with 16 4 sites (o). Statistical 
errors are comparable to the size of the dots. The continuous line represents a best fit of the form a + bx c . Finite 
size scaling predicts that all points should lie on the same universal curve. At k c — 0.0637 the scaling plot gives the 
value v = 0.333. 



depend on any fitting procedure, or at least analyzes a different and complementary set of data. By studying 
the dependence of averages on the physical size of the system, one can independently estimate the critical 
exponents. 

Fig. 4. shows a graph of the scaled curvature IZ(k) L A ~ X I V for different values of L = 4, 8, 16, versus the 



scaled coupling (k c — k~)L x l v . If scaling involving k and L holds according to Eq. ( 2.19 ), with xo = 1 — 4^ the 



scaling dimension for the curvature, then all points should lie on the same universal curve. From Eq. (2.19) 
with t ~ k c — k and xq = —5 = 1 — 4z/, one has 



n(k,L) = l-^- x i v [ n ((jfe c - fc) + 



(8.81) 



where uj > is a correction-to-scaling exponent. The data supports well such scaling behavior, and provides 
a further stringent test on the value for v, which appears to be consistent, within errors, with 1/3. 

Fig. 5. shows explicitly the size dependence of the average curvature. For small k the volume dependence 
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Fig. 5. Volume dependence of the average curvature, for L = 4, 8, 16, and (from top to bottom) k=0. 040, 0.045, 0.050 
and 0.055. 



is small, and gradually increases towards the critical point. Such a trend is in agreement with the expectation 
that the correlation length £ is growing as one approaches the critical point, leading to a more marked volume 
dependence. For fixed k ^ k c one expects on the four-torus 

n L (k) ~ Kooik) + Am(fc) 1 / 2 L~ 3 / 2 e - m{k)L + ••• , (8.82) 

L>l/m(fe) 

where L — V 1 ^ 4 is the linear size of the system and m = is the lightest mass in the theory. Combining 
and averaging the estimates from correlations Q, potential jL2| and finite size corrections to the average 
curvature one can in fact estimate the magnitude of this mass directly. One obtains m ~ 0.81 (k c — fc) 1 / 3 , 
giving a correlation length of about two lattice spacings at k = 0.050. 

The value of k c itself should depend on the size of the system. Indeed such a dependence is found when 
comparing k c (as obtained from the algebraic singularity fits discussed previously) on different lattice sizes. 
One writes 

k c {L) - fc c (oo) + c L~ 1/v + ■■■ . (8.83) 

L — >oo 

Fig. 6. shows the size dependence of critical coupling k c as obtained on different size lattices. In all three 
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Fig. 6. Volume dependence of the critical coupling k c , as determined from the singularity in the average curvature, for 
lattices with L — 4,8, 16. The points labelled by o are obtained assuming v — 1/3; in both cases the lines represent 
simple fits of the type k c = a + b/L 3 . 



cases k c (L) is first obtained from a fit to the average curvature of the form lZ(k) = A (k c — k) s for k > 0.02. 
Furthermore, if one assumes v — 1/3 and extracts k c from a linear fit to 1Z 3 , then the variations in k c for 
different size lattices are substantially reduced (points labelled by circles in Fig. 6.). Due to the few values 
of L it is not possible at this point to extract an estimate for v from this particular set of data. But since v 



is close to 1/3, it makes sense to use this value in Eq. ( p. 83 ) at least as a first approximation. 

Fig. 7. shows a plot of the average curvature TZ(k) versus reduced coupling k c — k, for several values of 



a, the higher derivative coupling of Eq. ( 4.26 ). a = corresponds to the pure Regge action with no explicit 
higher derivative lattice contribution, for which the path integral is still well defined (at least for sufficiently 
small \k\), since the deficit angles are bounded, and the edge lengths fluctuate around some average value, 
which is determined by the interplay of the measure and the cosmological constant term. Alternatively, one 
can think of the fluctuations in the conformal mode as becoming bounded (again at least for sufficiently 



small | A: |) when a momentum cutoff of order ir/V<l 2 > is dynamically generated. 

The slope of each straight lines determines the critical exponent 5 = Ay — 1, and it seems clear from the 
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Fig. 7. Average curvature 1Z versus reduced coupling fc c — k, on a log-log scale. From top to bottom, a — 
0,0.0005,0.005,0.02,0.1, with a the higher derivative coupling. Statistical errobars are comparable to the size of 
the dots. The slope of each straight lines determines the critical exponent 8 = 4u — 1. The slope is noticeably smaller 
for a = 0, suggesting that the higher derivative terms mask the true critical behavior up to very small fc c — k. 



graph that the slope is noticeably smaller for a = 0, suggesting that the higher derivative terms mask the 
true critical behavior up to very small k c — k (it was already noted in that for a = the assumption of 
an algebraic singularity for the average curvature leads to a value for the curvature exponent which is much 
smaller than the estimate for a > 0, namely S w 0.30(4)). 

Indeed it seems that one of the effects of the higher derivative terms is to push the region of instability 
towards smaller and smaller values of k c — k, until it becomes numerically undetectable. But we would argue 
that it is only close to this region that the correct continuum behavior is recovered. The situation is similar 
to what happens in the weak field expansion and perturbation theory: higher derivative terms do not cure 
the instability problems in the physically relevant region of small momenta and large correlation lengths. 

Fig. 8. shows a plot of the curvature amplitude An versus the higher derivative coupling a. The rapid 
growth close to a = is consistent with an expected catastrophic instability for a < (wrong sign for higher 
derivative terms). 
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Fig. 8. Curvature amplitude An versus the higher derivative coupling a. The amplitude increases rapidly as a 
approaches zero, the pure Einstein- Regge limit. 
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Fig. 9. Critical exponent v computed from the average curvature, versus the higher derivative coupling a. Note the 
small errorbar on the recent value for v at a = 0. For a > the higher derivative terms tend to mask the true critical 
behavior, which requires k c — k -C a~ . 



A compilation of previous estimates for is, together with the new value at a = 0, is shown in Fig. 9. 
There seems to be a clear trend toward smaller values as a approaches zero, the Einstein-Regge limit. While 
the Einstein action contribution becomes the dominant one at large distances, this is no longer the case at 
intermediate distances in the presence of the higher derivative terms. One concludes that for a > the 
higher derivative terms tend to mask the true critical behavior, which requires k c — k <C a . 



Fig. 10. shows a plot of the average volume per site < V>, in units of the average edge length y/<l 2 >- 
The curve is a fit of the form a + b(k c — k) c , and suggests a rather sudden drop of the average volume in 
the vicinity of the critical point. A non-analiticity in < V> at k c is in fact consistent with the sum rule of 



Eq. (5.45), which suggest that the singular behavior in the average curvature lZ(k) and the average volume 
<V> (k) are simply related. Typically, the sum rule of Eq. ( 5.45 ) is satisfied to one part in 10 3 or better. 



As can be seen from Fig. 10, close to the transition at k c the average volume per site expressed in 
units of the average lattice spacing, <V> / < I 2 > 2 shows only a weak singularity when the critical point 
is approached from the smooth phase (fc < k c ), and tends to a finite value. On the other hand, in the 
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rough phase (k > k c ) the volume per site seems to approach smaller and smaller values as the lengths of 
the runs are extended. In fact it would seem that in the rough phase the volume per site can be made to 
approach zero, at least for some simplices. One refers therefore alternatively to this phase as the collapsed 
or polymer-like phase, since its effective dimension is two i. Furthermore the relaxation times in the rough 
phase become exceedingly long, with the system getting stuck in some degenerate, spike-like configurations 
without being able to get out of it again. 

It seems difficult to see how the collapse of the simplices could be averted by choosing a different lattice 
structure (for example a random lattice), since its properties seem to be unaffected by changes in the 
measure or the action, at least to the extent they have been investigated. Indeed the collapsed, polymer-like 
phase appears even in the simplest models based on a regular tessellation of the four-sphere |Q, || . From a 
continuum point of view, the existence of such a pathological phase is not unexpected, and is interpreted 
as reflection of the unbounded fluctuations in the conformal mode expected for sufficiently large k. Indeed 
unbounded fluctuations in the conformal mode in the continuum correspond to rapid fluctuations in the 
simplicial volumes, and this is precisely what is observed on the lattice for k > k c , namely a rapid variation 
of simplicial volumes when going from one simplex to a neighboring one. 



Fig. 11. shows a plot of the average edge length \/<l 2 >. The curve is a fit of the form a + b(k c — k) c , 
and suggests a rapid increase in this quantity towards the critical point at k c . Indeed as the critical point 
is approached the number of fairly small and fairly large edge lengths proliferate, leading to an increasingly 
wide edge length distribution. 



Fig. 12. shows a plot of the average curvature fluctuation Xii(k) defined in Eq. (4.31). At the critical 
point the curvature fluctuation diverges, by definition. As in the case of the average curvature lZ(k) analyzed 
previously, one can extract the critical exponent S and k c by fitting the computed values for the curvature 
fluctuation to the form (see Eq. ( 6.56D ) 

Xn(k) ^ A XK (k c - fc)~ (1_a) . (8.84) 

As for the curvature itself, it would seem unreasonable to expect that the computed values for 1Z are 
accurately described by this function even for small k. Instead the data is fitted to the above functional form 
for either k > 0.02 or k > 0.03 and the difference in the fit parameters can be used as one more measure for 
the error. Additionally, one can include a subleading correction 

Xn(k) ^ -A xn [k c ~k + B(k c - k) 2 ] - {1 ~ 5) , (8.85) 



3 An elementary argument can be given to explain the fact that the collapsed phase for k > k c has an effective dimension of 
two, as was found in |12j. The instability is driven by the Euclidean Einstein term in the action, and in particular its unbounded 
conformal mode contribution. As the manifold during collapse reaches an effective dimension of two this term turns into a 
topological invariant, unable to drive the instability further to a still lower dimension. 
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k = 1/8ttG 

Fig. 12. Curvature fluctuation on lattices with 4 4 (□), 8 4 (A) and 16 4 (o) sites. The thin-dotted, dotted and 
continuous lines represent best fits of the form Xiz(k) = A (k c — fc)~' 1 ~' 5 ' ) for k > 0.02. 

and use the results to further constraint the errors on A XK , k c and 5 = 4v — 1. 

The values for S and k c obtained in this fashion are consistent with the ones obtained from the average 
curvature TZ(k), but with somewhat larger errors, since fluctuations are more difficult to compute accurately 
than local averages, and require much higher statistics. Using these procedures one obtains on the lattice 
with 16 4 sites 

k c = 0.0636(30) v = 0.317(38) . (8.86) 

Fig. 13. shows a graph of the inverse curvature fluctuation Xn{k) on the 16 4 -site lattice, raised to power 
3/2. One would expect to get a straight line close to the critical point if the exponent for xn{k) is exactly 
—2/3. The numerical data indeed supports this assumption. The computed data is consistent with linear 
behavior for small k > 0.02, providing further support for the assumption of an algebraic singularity for 
Xn{k) itself, with exponent close to —2/3. Using this procedure one finds on the 16 4 -site lattice 

k c = 0.0641(17) , (8.87) 

which is completely consistent with the value obtained from 1Z 3 (see Fig. 3. and related discussion), and 
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Fig. 13. Inverse curvature fluctuation raised to the power 3/2, on the 16 4 (o) lattice; data is scaled by a factor of 
xlOO. The straight line represents a linear fit of the form A (k c — k). The location of the critical point in k is 
consistent with the estimate obtained from the average curvature, but with a somewhat larger error. 



suggests again that the exponent v must be close to 1/3. 

Fig. 14. shows the results for the logarithmic derivative of the average curvature IZ(k), obtained from 
the data shown in Figs. 3. and 12. From the definition of the average curvature 1Z and curvature fluctuation 
(Eqs. ( 4.29| ) and ( 4.31 )), and the fact that both are proportional to derivatives of the free energy F with 



respect to k (Eqs. ( 4.34 ) and ( 4.35 )), one notices for the ratio 



—}Z(kj { dk lnZL)/i dP hjZL) ~ ~dk \ ~dk 



(8.88) 



The assumption of an algebraic singularity in k for TZ and xiz (Eqs. (6.55) and (6.56)) then implies that the 
logarithmic derivative as defined above has a simple pole at k c , with residue 5 — Au — 1 

*J£> Xnjk) „ _J_ 

K{k) k->k c k-k c ' y ' ' 

with the critical amplitude dropping out of this particular expression. The above result is general and does 
not rely on k being real. This suggests that in principle the method of Pade rational approximants (which 
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Fig. 14. Inverse of the logarithmic derivative of the average curvature lZ(k). The straight line represents a best fit of 
the form A (k c — k) for k > 0.02. The location of the critical point in k is consistent with the estimate coming from 
the average curvature TZ. From the slope of the line one computes directly the exponent v. 

applies only to meromorphic functions) can be employed to locate singularities in xti(k), even for complex 
k [[l3| [ll) . Using this method on the 16 4 lattice one finds 



0.0635(11) v = 0.339(9) 



.90) 



It is encouraging that the above estimates are in good agreement with the values obtained previously using 
the other methods. 

Fig. 15. shows a graph of the scaled curvature fluctuation xtz (k) / L 2 / u ~ 4 for different values of L = 4, 8, 16, 
versus the scaled coupling (k c — k)L x l v . If scaling involving k and L holds according to Eq. ( J2.19| ), with 
t ~ k c — k and xq = 1 — 5 = 2 — 4v then all points should lie on the same universal curve. From the general 



Eq. (2.19) one expects in this particular case 

Xn(KL) = L 2 ^- 4 [ X n{{k c -k)L^) + O(L^) 



(8.91) 



where u) > is again the correction-to-scaling exponent. Again the data supports such scaling behavior, and 
provides a further estimate on the value for v, close to 1/3. 
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Fig. 15. Finite size scaling behavior of the scaled curvature fluctuation versus the scaled coupling. Here L — 4 for 
the lattice with 4 4 sites (□), L = 8 for the lattice with 8 4 sites (A), and L = 16 for the lattice with 16 4 sites (o). 
The continuous line represents a best fit of the form l/(a + bx c ). Finite size scaling predicts that all points should 
lie on the same universal curve. At k c = 0.0637 the scaling plot gives the value v — 0.318. 
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Fig. 16. shows a plot of the curvature fluctuation xiz versus the curvature 1Z. If the curvature approaches 
zero at the critical point where the curvature fluctuation diverges, one would expect the curvature fluctuation 
to diverge at TZ = 0. One has 

X n(K) fc ~ A\TZ\ (1 - S)/S ~ A \ n \(4»-2)/(4v-i) (8 92) 

An advantage of this particular combination is that it does not require the knowledge of k c in order to 
estimate v. Using all points corresponding to k > 0.02 one finds 

v = 0.328(6) . (8.93) 

The error on v can be estimated, for example, by using a more elaborate fit of the type 

Xn n ~ Q A | K + BK 2 |( 4 "- 2 )/( 4 — i) . (8.94) 

For v = 1/3 the exponent becomes equal to —2, and one has the simple result 

X n n - Q A \K\- 2 . (8.95) 

One concludes that the evidence supports a vanishing curvature at the critical point, where the curvature 
fluctuation xiz and the correlation length £ diverge. This result is further supported by the consistency of 
the values for k c obtained independently from IZ(k) and XTt(k) (Figs. 2,3,4,12,13,14 and 15). 

As an independent measure of the fluctuation one can also investigate the behavior of the edge length 
fluctuation defined as 



1 ( Ni Ni 

< k ) - wA <E<?) a >- 

1 l i=l f=l 



\ ,-•(/.■) ^ — i >■ <}jjr •:. > J:>^ (k c ~k)^ , (8.96) 

where 7 is a critical exponent. Using an analysis similar to what is done for the curvature and curvature 
fluctuation, on the 16 4 lattice it is found to diverge at 

k c = 0.0609(23) (8.97) 

in agreement within errors with the previous values quoted for k c . One would expect such a fluctuation 



to be related to the fluctuations in the local volumes, and, by the sum rule of Eq. ( 5.45 ) which relates the 
fluctuations in the volume to fluctuations in the curvature, one would expect 7 = 1 — 5 = 2 — Av. The 
numerical results for gamma have larger errors but give values between 0.46 and 0.85, certainly consistent 
with a value of 7 = 2/3 for v = 1/3. 

Finally Fig. 17. summarizes the known information about the phase diagram in the k-a plane. The 
continuous line separates the smooth phase with small negative curvature from the rough, polymer-like 
phase. 
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Fig. 17. Phase diagram for the model in the k — a plane. A critical line separates the smooth, strong coupling, phase 
from the rough, weak coupling, phase. The dotted line denotes the pure Einstein theory, without higher derivative 
terms. 
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Method 


k c 


V 


TZ vs. k 


0.0630(11) 


0.330(6) 


TZ A vs. k 


0.0639(10) 


- 


Xn vs. k 


0.0636(30) 


0.317(38) 


3/2 , 

Xn vs. k 


0.0641(17) 


- 


Xn/mn) vs. k 


0.0635(11) 


0.339(9) 


Xn vs. TZ 




0.328(6) 


Xi 2 vs. fe 


0.0609(23) 


7 = 0.46(8) 


X/a vs. TZ 




7 = 0.54(7) 


1Z FS scaling 




0.333(2) 


Xn FS scaling 




0.318(10) 


FS scaling 




7 = 0.85(6) 



Table I: Summary table for the critical point k c and the critical exponent v, as obtained from the largest lattice with 
16 4 sites. The last three entries assume a critical point at k c — 0.0636. 

Table I summarizes the results obtained for the critical point k c — 1/8ttG c and the critical exponent v. 
From the best data (with the smallest statistical uncertainties and the least systematic effects) one concludes 

k c = 0.0636(11) v = 0.335(9) , (8.98) 

which suggests v = 1/3 for pure gravity @. 



9 Critical Exponents and Phenomenoloy 

In this section some consequences of the results presented above will be discussed, with ultimately an 
eye towards possible physical applications. Naively one would expect simply on the basis of dimensional 
arguments that the curvature scale gets determined by the correlation length 

TZ - IIS 2 , (9.99) 
n^o 

but one cannot in general exclude the appearance of some non-trivial exponent. 

In the previous section arguments have been given in support of the value v = 1/3 for pure gravity. From 



Eqs. ( 6.58 ) relating the average curvature to the correlation length one has 



^(0 ~ 9 wXTT ) (9- 10 °) 



4 The value v = 1/3 does not correspond to any known field theory or s tatistical mechanics model in four dimensions. For 
dilute branched polymers it is known that v = 1/2 in three dimensions Ma], and v = 1/4 at the upper critical dimension d = 8 
|47|, so one would expect a value close to 1/3 somewhere in between. I thank John Cardy for a discussion on this point. 
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and the correct dimension for the average curvature 1Z have been restored by supplying appropriate powers 
of the ultraviolet cutoff, the Planck length / p = y/G. One notices that close to two dimensions the exponent 
of £ indeed approaches 2, since v ~ l/(d— 2), and the classical result is recovered. 
For v = 1/3 in four dimensions i one then obtains the remarkably simple result 

n ® k ~ k .ih ■ (9 - 101) 

An equivalent form can be given in terms of the curvature scale Hq, defined through R = —12Hq, and which 
has dimensions of a mass squared. One has close to the critical point 

Hi = C H Mp m , (9.102) 

where /ip = \/\/G is the Planck mass, m = l/£ is the inverse gravitational correlation length, and Ch ~ 4.9 
a numerical constant of order one; the value for C is extracted from the known numerical values for 1Z and 
m close to the critical point at k c . 

One can raise the legitimate concern of how these results are changed by quantum fluctuations of matter 
fields. In the presence of matter fields coupled to gravity (scalars, fermions, vector bosons, spin-3/2 fields etc.) 
one expects the value for v to change due to vacuum polarization loops containing these fields. A number of 
arguments can be given though for why these effects should not be too dramatic, unless the number of light 
matter fields is very large. First, in the case of a single light scalar field the vacuum polarization effects are 
so small that they are barely detectable in the numerical evaluations of the path integral 0] . Furthermore 
one notices that to leading order in the 2 + e expansion the exponent v only depends on the dimensionality 
of space-time, irrespective of the number of matter fields and of their type ]4lf| . Finally one can compute for 
example the effects of scalar matter fields on the one-loop beta function in the 2 + e expansion for gravity, 



and finds /3q — (2/3) (25 — n/) where n/ is the number of massless scalar fields 41 . Thus unless n/ is large, 
the matter contribution is quite small even to next-to-leading order in the 2 + e expansion. The present 
evidence would therefore suggest that the approximation in which vacuum polarization effects of light matter 
fields are neglected should not be too unreasonable. 

It seems natural to identify Hq with either some (negative) average spatial curvature, or possibly with 
the Hubble constant determining the macroscopic expansion rate of the present universe [ fl2[ |l5f| . In the 
Friedmann- Robertson- Walker model of standard cosmology |5Cj on has for the Ricci scalar 



5 For all scalar field theories (spin s = 0) in four dimensions it is known that v = 1/2, while for the compact Abelian U(l) 
gauge theory (s = 1) one has v = 2/5 [MSI. The value v = 1/3 for pure gravitation (s = 2) in four dimensions is then consistent 
with the simple formula v = 1/(2 + s/2). 
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where R(t) is the FRW scale factor, and k — 0, ±1 for spatially flat, open or closed universes respectively. 
Today the Hubble constant is given by H$ = (R/R)f Q , but it is eventually expected to show some slow 
variation in time, and its characteristic length scale cHq 1 10 28 cto today is comparable to the present 



extent of the visible universe. Under such circumstances from Eq. (9.102) one would expect the gravitational 
correlation £ to be significantly larger than cH^ 1 . A potential problem arises though in trying to establish 
a relationship between quantities which are truly constants (such as the ones appearing in Eq. fl9J02T) ), and 
Hq which most likely depends on time. 1 In any case it is clear that some of these considerations are in fact 
quite general, to the extent that they rely on general principles of the renormalization group and are not tied 
to any particular value of is, although v = 1/3 clearly has some aesthetic appeal. Additional cosmological 
and astrophysical arguments and proposed tests can be found in a recent paper |5lJ . 

One further observation can be made regarding the running of G. Assuming the existence of an ultraviolet 



fixed point, the effective gravitational coupling is given by Eq. (7.66) for "short distances" r <C £, but now 
with an exponent v = 1/3, 

G(r) = G(0) [1 + c(r/0 3 + 0((r/£) 6 )] , (9.104) 
with c a calculable numerical constant of order one. The appearance of £ in this equation, which is a very 



large quantity by Eq. (9.102), suggests that the leading scale-dependent correction, which gradually increases 
the strength of the effective gravitational interaction as one goes to larger and larger length scales, should 
be extremely small. @ 

It is only for distances comparable to or larger than £ that the gravitational potential should start to 
weaken and fall off exponentially, with a range given by the gravitational correlation length £, 

V(r) ~ - G(r) MlM2 . (9.105) 

r >f r 

In many ways these results appear qualitatively consistent with the expected behavior of the tree-level 
graviton propagator in anti-de Sitter space J52|, [5j| . In the real world the range £ must be of course very 
large. From the fact that super-clusters of galaxies apparently do form, one can easily set an observational 
lower limit £ > 10 25 cm. 

It is unclear to what extent gravitational correlations can be measured directly. From the definition of 



the curvature correlation function in Eq. (6.53) one has for "short distances" r <C £ and for the specific value 
v = 1/3 the remarkably simple result 

<^R(x) ^gR(y)8(\x-y\-d)> c ^ ^ , (9.106) 



6 The only exception being the steady state cosmological models, where H is truly a constant of nature. These models are 
not favored by present observations, including detailed features of the cosmic background radiation. 

7 And suggests that the deviations from classical general relativistic behavior for most physical quantities is in the end 
practically negligible. 
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with A a calculable numerical constant of order one. One can contrast this behavior with the semiclassical 
result attained close to two dimensions (and which incidentally coincides with the lowest order weak field 
expansion result |3q| ), which gives instead for the power the value 2(d— 1/v) ~ 2(d— (d— 2)) ~ 4, as expected 
on the basis of naive dimensional arguments (R ~ d 2 h). 

If one considers the curvature R averaged over a spherical volume V r = 47rr 3 /3, 

VgR = [ d 3 xy/g(x,t)R(x,t) (9.107) 

*V JV r 

one can compute the corresponding variance in the curvature 

[S(V9 R)f = ^ J d3£ J d3 v < R@) Va W) >c = |^ ■ (Q.ios) 

As a result the r.m.s. fluctuation of y/gR averaged over a spherical region of size r is given by 

2 r 

while the Fourier transform power spectrum at small k is 

4-TT 2 A 1 

p% = lv«l 2 = ~ k ■ ( 9 - n °) 

One can use Einstein's equations to relate the local curvature to the (primordial) mass density. From 
Einstein's field equations 

Rp.v - \g^R = 8wT„ v (9.111) 

for a perfect fluid 

Tpv = pg^u + (p + p^Uu (9.112) 
one obtains for the Ricci scalar, in the limit of negligible pressure, 

R(x) w 8nG p(x) . (9.113) 

As a result one expects for the density fluctuations a power law decay of the form 

< p(x) p(y) > c ~ . (9.114) 

k-vK£ \x~y\ 2 

Similar density correlations have been estimated from observational data by analyzing known galaxy number 
density distributions, giving a value for the exponent of about 1.77±0.04 for distances in the 10/cpc to lOMpc 
range Q. 
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10 Concluding Remarks 



Numerical simulation methods combined with modern renormalization group arguments and finite size 
scaling can provide detailed information on non-perturbative aspects of a lattice model of quantum gravity. 
It has been known for some time that the lattice model has two phases, only one of which is physically 
acceptable. In this work we have described in some detail the properties of the latter smooth phase, and 
provided quantitative estimates for the critical point, the scaling dimensions and the behavior of correlations 
at distances large compared to the cutoff. In spite of the fact that the Euclidean theory becomes unstable 
as one approaches the critical point at k c , it is still possible to determine by a straightforward analytic 
continuation the physical properties of the model in the vicinity of the true fixed point, defined as the point 
where a non-analiticity develops in the strong coupling branch of Z^k), and where scaling implies that the 
physical correlation £ diverges. 

If this prescription is followed, an estimate for the non-perturbative Callan-Symanzik beta function in 
the vicinity of the fixed point can be obtained, to leading order in the deviation of the bare coupling from its 
critical value. The resulting scale evolution for the gravitational constant is then quantitatively quite small, 
if one assumes that the scaling violation parameter is related to an average curvature and its characteristic 
scale H - Its infrared growth, consistent with the general idea that gravitational vacuum polarization effects 
cannot exert any screening, suggests that low energy properties of quantum gravity are inaccessible by 
weak coupling perturbation theory: low energy quantum gravity is a strongly coupled theory. On a more 
quantitative side, as pointed out in the discussion there are a number of attractive features to the pure 
gravity result v = 1/3, including a simple form for the curvature correlations at short distances. 

It seems legitimate to ask the question whether the present lattice model for quantum gravity provides 
any insight into the problem of the cosmological constant. The answer is both yes and no. To the extent 
that a naive prediction of quantum gravity is that the curvature scale should be of the same order of the 
Planck length, 1Z ~ 1/G, the answer is definitely yes. Indeed it can be regarded as a non-trivial result of the 
lattice models for gravity that a region in coupling constant space can be found where space-time is stiff and 
the curvature can be made much smaller than 1/G. In fact the evidence indicates that the average curvature 
1Z vanishes at the critical point k c . And this is achieved with a bare cosmological constant A which is of 
order one in units of the cutoff. Phrased differently, the dimensionless ratio between the renormalized and 
the bare cosmological constant becomes arbitrarily small towards the critical point. 

At the same time the effective long distance cosmological constant is non- vanishing and of order l/£, 
and the value zero is only obtained when £ is exactly zero, which happens only at the critical point k c . 
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Thus to make the effective cosmological constant small requires a fine tuning, in the sense that the bare 
coupling k c — k has to be small. But since the correlation length determines the corrections to the Newtonian 
potential (and in particular its eventual decrease for large enough distances) , it would seem unnatural to have 
a short correlation length £: in such a world there would be no long-range gravitational forces, and separate 
space-time domains would have decoupled fluctuations. From this perspective, long range forces and a small 
cosmological constant go hand in hand. Quantum fluctuation effects show that hyperbolic space-times with 
small curvature radii cannot sustain long-range gravitational forces, at least in this model. 
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